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Abstract: 

The production of top quark pairs accompanied by a hard gluon in e + e~ annihilation is 
studied including next-to-leading order corrections in the strong coupling. At leading order, 
the fraction r of tig events with respect to all ti events is computed analytically as a function 
of the minimal gluon energy. Next-to-leading order results for r are given for center-of-mass 
energies of 0.5 and 1 TeV. We further calculate the differential distribution of r with respect 
to several variables, including the top quark energy and the ti invariant mass. We then 
investigate how our results depend on the choice of the renormalization scheme for the top 
quark mass by comparing results expressed in terms of either the pole mass or the MS mass. 
Finally we estimate the sensitivity of the fraction r on the value of the running top quark 
mass at a scale of 1 TeV. 
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1 Introduction 



A future high-energy e + e~ collider will open up a new domain for the experimental in- 
vestigation of the fundamental interactions between elementary particles. In particular, 
the production of top quark pairs in a clean environment will offer many interesting op- 
portunities both to perform precision tests of Quantum Chromodynamics (QCD) and to 
probe possible deviations from the standard theory. For example, the total cross section 
for e + e~ — > ti has been computed recently to next-to-next-to-leading order in the strong 
coupling a s both in the threshold region |TJ, |3|] and also far above threshold M. The 
results obtained for the threshold region will allow for a very precise determination of the 
top quark mass [§, [| 0] • 

The subject of this paper is the production of top quark pairs above threshold accom- 
panied by (at least) one additional parton with a hard momentum. To be more explicit, we 
are interested in the reaction 



e + (p + ) + e"(p_) - t(k t ) + + X(k x ), (1) 

and we require that the parton(s) X carry a large momentum kx- At order a s , the additional 
parton can only be a gluon, and at order a 2 s we have the possibilities X = g, gg, qq. By 
determining the properties of a final state with three or more partons one can test the 
strong interactions in great detail as is well known from the extensive studies of jets in 
e + e~ annihilation. For example, by comparing the cross section for (H) (which will be 
precisely defined in the next section) to the inclusive ti production cross section, one could 
measure a s . Any deviation from the standard determination of a s using event samples 
containing light quark jets would indicate a violation of the "flavour independence" of the 
strong interactions — i.e., would point towards new physics phenomena. A specific example 
of non-standard interactions that could be probed by reaction ([I]) are possible anomalous 
couplings of the top quark to photons, Z-bosons and gluons. In || it has been shown that a 
large anomalous chromomagnetic tig coupling would modify the gluon energy spectrum in 
e + e~ — > tig. Furthermore, symmetry tests can be performed utilizing the richer kinematic 
structure of the final state in ([]]). These contain tests of the CP symmetry |§ and the search 



for final state rescattering effects using triple momentum correlations [HI. Both the search 



for heavy quark anomalous couplings |TT| and the symmetry tests [0 have been shown to 
be experimentally feasible in the case of bbg production at the Z resonance, and it will be 
interesting to see whether similar studies are possible with top quarks. 

For any of the above studies it is mandatory to analyse reaction ([!]) at next-to-leading 
order in a s , which is the topic of this paper. It is organized as follows: We start in section |2| 
by studying reaction (|TJ) at leading order. We derive an analytic formula for the fraction r of 
tig events with respect to all ti events as a function of the minimal gluon energy. In section |3| 
we discuss the QCD corrections to reaction ([[]). We evaluate the fraction r at next-to- leading 
order for center-of-mass energies of 0.5 and 1 TeV. We further compute the differential 
distribution of the cross section for reaction (|l|), normalized to the total ti cross section, with 
respect to the top quark energy, the ti invariant mass, and the cosine of the angle between 



1 



the t and i. In section |] we study the dependence of our results on both the renormalization 
scale and the renormalization scheme employed to define the top quark mass. We first 
review the dependence of the total ti cross section on the mass renormalization scheme by 
expressing the result to order a s in terms of either the perturbative pole mass or the MS 
(running) mass. We then study in the same fashion the scheme and scale dependence of the 
fraction r. Finally we estimate, for a cm. energy of 1 TeV, the sensitivity of the fraction r 
on the value of the running top quark mass. 



2 Analysis at leading order 

In this section we discuss in some detail reaction flU) at leading order (LO) in at s , i.e., the 
production of a ti pair together with a single gluon, 

e+(p+) + e~(p_) -> t{kt) + t{k) + g{k g ). (2) 

We start by defining the following dimensionless variables: 
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%tt — ~ — 1 ~~ Xg, (3) 

where y/s is the center-of-mass energy, k = p + + p_, and E t f y9 are the energies of the 
final state particles in the cm. system. The cross section for reaction (|2]) develops a soft 
singularity as E g — >• 0. An infrared finite cross section may be defined by demanding 

•I'd •'-cut ; (4) 

where x cnt is some preset number. This condition avoids the region of phase space where the 
gluon becomes soft. Due to the finite mass of the top quark no collinear (mass) singularities 
arise and thus the above condition is sufficient to define an infrared finite cross section. The 
requirement (||) does not lead to a finite cross section for e + e~ — > qqg, where q is a massless 
quark, because the latter cross section is also singular when the gluon is collinear to the 
( ant i-) quark. 

We want to study the fraction r of tig events for which x g > x cn t with respect to all ti 
events. Since x g is no useful variable for final states with four or more partons (which will 
be relevant at higher orders in a s ), we use instead the scaled ti invariant mass square x t i 
defined in (|3|), i.e. replace the condition x g > x cut by the (at LO equivalent) condition: 

1 — Xti > ^cut- (5) 
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The fraction r is defined as 



/ , _cr (e+e -> ttX; 1 - x it - > x cut ) _ a 3 (x cut ) , . 

o- tot (e+e- -> tt) o- tot 

At LO we have X = g and write: 

r(x cut ) = ss fM* cut ) + O (a?) • (7) 

Here, a s /(27r)cx3(a; cu t) denotes the LO cross section for e + e~ — ► tig with 1 — x t t > x cut and 
of ot is the LO inclusive cross section for e + e~ — > ti. The coefficient A(x cut ) may be written 
as follows: 



A(x. 



cut i 



2N C CF^pt f dx f dxF 1 (x,x)Q(l - x t t - x cut )Q(l -cos 2 9 t t) 

a tot JO JO 



where Cp = (Nq — l)/(2Nc) and Nc = 3 is the number of colours. The point cross section 
apt reads 

apt = cr(e e ^7 -> fi fi ) = , (9) 

and the LO total cross section of ot can be expressed in terms of the dimensionless mass 
variable, 



s 



(10) 



as follows: 

a t ot = N c a pt VT^ [c v (l + 2z) + c A (l - Az)] . (11) 

The coupling factors cy,A appearing in (|ll|) are given explicitly by: 

Cv = Q 2 t r +2g v Q t Rex(s) f z + g t2 \ X (s)\ 2 f\ 

CA = 9 t A 2 \x(s)\ 2 f zz , (12) 

with 

H = i- A_A + , 

f™ = (1 _ A_A+)(^ 2 + g\ 2 ) - 2(A_ - X+)g v g A , 

fiz = _(i _ A_A+)<£ + (A_ - X + )g% (13) 

Here, Q t = 2/3 is the electric charge of the top quark, and g A v denote the axial and vector 
couplings of the fermion /. In particular, g v = — | + 2sin 2 ^ H /, g\ = —\ for an electron, 



and 9y — 2~ fsin 2 ?^, q\ = \ 101 a top quark, where $w is the weak mixing angle. The 
function x{ s ) reads 



4 sin ^vy cos 2 $w s ~ m z + ^ m z^z ' 



(14) 



where mz and stand for the mass and the width of the Z boson. Finally, A_ (A+) denotes 
the longitudinal polarization of the electron (positron) beam. 
The function F\ may be decomposed as: 



Fi(x,x) = c v F^(x,x) + c A F 1 (x,x). 



(15) 



Using the abbreviation 



B 



(l-x)(l-x) 



(16) 



we have |13 

F%{x,x) = b{ 



x 2 + x 2 



[2x g + ((1 - x) 2 + (1 - x) 2 ) B] - 2z 2 x 2 g By (17) 



Ff(x, x) = F^ v (x, x) + zB^(x + x) 2 - 10(1 - x g ) + 6za£s}. (18) 

The cosine of the angle between the top quark and antiquark appearing in (|8]) can be 
expressed in terms of the top quark mass and the scaled energy variables: 

xx — 2(1 — x a ) + Az , . 

cos Bit = = (19) 

a/(x 2 - Az) (x 2 - Az) 

The Heaviside function 0(1 — cos 2 9 t i) defines the kinematically allowed region in the (x, x) 
plane when no cuts are applied. For m t = 175 GeV and the two cm. energies y^i = 0.5 
TeV and y/s = 1 TeV this region is depicted in Fig. [IJ The cuts on the scaled gluon energy 
a^cut = 0.1, 0.2 are indicated as lines, and the enveloping triangle is the kinematically allowed 
region for massless quarks. With a cut x cut on the gluon energy, the kinematically allowed 
area P(x cut , z) is: 



P(x cut , z)= dx dx 0(1 - x t t ~ x cut )0(l - cos 2 Q t t) 
Jo Jo 



'0 

= 0(1 -Az- x cut ] 

where 



2z(l - z) In (p) + 1(1 - s cut )(l + x cut + 2z)- 



1 - >/l - 4z/(l - x cnt ) 

l + y/l- Az/{1 - X cnt ) 



(20) 



(21) 
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When the kinematic boundary 1 — Az = x cnt is approached for a fixed value of x cnt , P{x C nt, z) 
vanishes like (1 — Az — x cut ) 3 ^ 2 . 

The function A(x cnt ) defined in (§) can be very easily computed numerically. This has 
the advantage that an implementation of additional and/or different cuts is straightforward. 
On the other hand, an analytic result for A(x cut ) is obviously desirable, and in fact it is 
possible to perform the twofold integral in (§) analytically. The result is: 

A(x cut ) = 2NcC F(Jvt Q{l -Az- x cut ) [c v A v (x cut ) + c A A A (x cut )} , (22) 

a tot 



where 



with 



A v (x cnt ) = --(1 + 4a; cut - x 2 cut + 8zx cut + 2z 2 ) ln(p) 

+ 7(1 -x cut )(ll -3x cut + 3Az)\^ - (1-Az 2 )g(u,p), 
A 1 + p 

A A (x cut ) = A v (x cut ) -z (1-Az- 12x cut - x 2 cut + 6z 2 ) ln(p) 

+ -(1 - a;^) (51 + x cut - 6z) - 6(1 - 2z)g(u,p)] , (23) 

Z 1 "T P 



g(uj, p) = ln 2 (u) - 2 In (^-^-) ( In( W ) + ^—^ 

\1 - PUJ J \ l+LU 2 

+ ln 2 (p) + 2 Li 2 (l - + 2 Li 2 (1-pu) (24) 



and 



to = . (25) 

l + ^/l^A^ 

Fig. y shows the fraction r defined in (|^) as a function of the cm. energy for m t = 175 
GeV, and three different values of x cut . Here and in all other numerical results of this paper 
we consider unpolarized electron and positron beams. The running of the strong coupling is 
taken into account in the curves, i.e. we use a s (p = y/s) with six active flavours. Our input 
value is a™ f ~ 5 (p = m z ) = 0.118, which is evolved up to p, — m t with five active flavours, 
and then converted into a™ f ~ 6 (m t ) using the so-called matching conditions of F° r 
example, we thus get a™ f ~ 6 (p = 500 GeV) = 0.0952. 

We close this section with a remark concerning the experimental distinction of ti events 
with a gluon radiated off the t or t from events in which the gluon is radiated off the b or b 



produced in the decays of the top quark pairs. It has been shown in |15] that the following 



lr To be more precise, the conversion is performed at the scale m t (mt), where m t is the running top quark 
mass. 
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two constraints efficiently select events where the gluon is produced in association with the 
top quark pair: 



By requiring that the invariant mass of the Wb system lies in the vicinity of the top quark 
mass, the probability that a highly energetic gluon jet (E g > T ( ^ 1.4 GeV) is emitted 
from the b or b is very small. Our computation may therefore be applied to describe events 
of the type e + e~ — > W + W~bbg that fulfil both conditions of (p6|) . Our "default" value for 
many of the results below will be x cut = 0.1, which corresponds to E g = 25 (50) GeV for a 
cm. energy yfs = 0.5 (1) TeV. 

3 Results at next-to-leading order 

At order a 2 s , we have to consider both virtual and real corrections to the process (||). The 
real corrections consist of the processes e + e _ — > ttgg, e + e~ — > tiqq (q = u, d, s, c, b). If \fs > 
4mt, the production of two ti pairs becomes possible. However, these rather spectacular 
events are extremely rare for the cm. energies considered below and contributions from the 
process e + e~ — > titi can therefore be neglected^]. 

As mentioned before, the condition 1 — x t i > x cut is equivalent, at LO, to requiring a 
minimal scaled gluon energy x g . For final states with four or more partons, an alternative 
condition is 2kkx/s > x cut , where in the cm. system 2kkx/s is equal to 2Ex/yfs, with 
Ex = y/s — E t — Ef. For a given cut x cut , this condition leads to larger contributions to 
K^cut) fro m the four-parton final states as compared to the definition employed in @, since 
2E X /Vs = 1 - x a + k\/s with k\ > 0. 

We renormalize the coupling in the modified minimal subtraction (MS) scheme. For the 
results of this section, the top quark mass is defined as the perturbative pole mass, i.e., 
the mass renormalization is carried out in the on-shell scheme. The result for the fraction 
r(x cut ) at a renormalization scale \x may be written to next-to- leading accuracy as 




(26) 



CtotO") 



2tt 



( 



2tc 



2 



K^cut; A 4 ) 



A(x cat ) + 



B(x cat ,fj) + 0(al) 



(27) 



with 




(28) 



2 Contributions from "secondary" t-quarks, originating from gluon splitting in e + e 
heavily suppressed. 



qqg, are also 
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^W = ^ + ^« + OW)- (29) 



The coefficient A(x cut ) is given explicitly in Eqs. ([22|), (p3|) , and the coefficient B(x cut ,fi) 
reads: 

p/_ ,,\ _ ojQEcut, AO gtotgfCgcut) /o n \ 

-D^cut^J — o / n \2 • 

^tot « ot ) 

The result for a\ ot is well known (see, e.g. [0]). For the computation of a\{x cnt ,^) we use 



the techniques and results of |17| and Ul3| . In these papers, the production of three jets 
involving b quarks was studied at NLO, taking into account the b quark mass. The matrix 
elements given there can be easily adapted to the case at hand. For technical details of the 
calculation, we refer the reader to [[13] . 



Figure 3 (4) shows the LO and NLO results for r as a function of x cut at = 0.5 TeV 
(y^i = 1 TeV). The renormalization scale is set to fi = y^. (The scale dependence of our 
results will be discussed in the next section.) In Table 1 we list the A and B coefficients for 
a sample of x cut -values. The numerical errors in the last digits of the B coefficients are also 
shownP]. For yfs = 0.5 TeV, the relative size of the QCD corrections, a s /(2%)B/A, varies 
between 56% (at a; cut = 0.02) and 36% (at x cut = 0.2). At s/s = 1 TeV, the QCD corrections 
are roughly constant as a; cut is varied and of the order of 30%. In Figs. 5, 6, and 7 we plot 
various distributions of the cross section a 3 (x cut , fi). All distributions are normalized to the 
total ti cross section <7 to t, and we set y/s = fi = 0.5 TeV, m t = 175 GeV, and x cut = 0.1. 
These values lead to the kinematic limits 0.7 < x < 0.9837, 0.1 < \—x t i < 0.51 for the scaled 
top quark energy x and the quantity l — x t t, respectively. (Recall that in leading order 1 — x t i 
is equal to the scaled gluon energy x g .) Fig. 5 shows the distribution l/a tot da 3 /dx. The 
distribution reaches its maximum around x ~ 0.925, which corresponds to E t « 231 GeV. 
We note in passing that the distribution with respect to the scaled top antiquark energy x 
is equal to the one shown due to charge conjugation invariance of the strong interactions. 
In Fig. 6 we plot the distribution l/a to t da 3 /d{l — x t i)- The size of the QCD corrections 
depends only weakly on 1 — x t t- (The seemingly large correction in the last bin is most 
probably a fluctuation due to the limited statistics of the Monte Carlo integration.) Fig. 7 
depicts the distribution with respect to the cosine of the angle between the top quark and 
antiquark. As one might have expected, this distribution is very sharply peaked close to 
cos# tt - = -1. 



4 Top quark pole mass versus running mass 

All results in the preceding section have been obtained by using the MS scheme for the 
coupling renormalization and by renormalizing the top quark mass in the on-shell scheme. 
In this section we study the dependence of our results on the mass renormalization scheme 
as well as on the choice of the renormalization scale \i. 

3 For example, 57.1(10) stands for 57.1 ± 1.0. 
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In order to understand the impact of a change of the mass definitions on our fixed order 
predictions, it is instructive first to consider the simple and well-known case of the inclusive 
rate a tot . We use the notation z on = (m° n ) 2 fs and z(ji) = m 2 (fi)/s, where m™ is the pole 
mass (for which we always use the value m™ = 175 GeV) and m t (/j,) is the MS mass at 
a scale \i. The latter is obtained from the renormalization group evolution with the input 
fn t (fi = frit) — 168 GeV. The relation between the two mass parameters is given by 

z on = z(fi) + Az(fi), (31) 

where, to order a s , 

z(ii) + O {a 2 s ) . (32) 
To order a s we thus havef]: 

= aUm) + A ^)^f W)) + ^oL&M) + O K) , (33) 

where da^ ot /dz can be obtained from ([TT]) . Fig. 8 (9) shows the total cross section in units 
of the point cross section @ as a function of the renormalization scale at y/s = 0.5 TeV 
(y/s = 1 TeV) for both mass definitions. When cr to t is expressed in terms of the pole mass, 
the dependence on /i is (to order a s ) only induced by the running of a s . Two important 
features are visible in the curves: First, within a large range of values for fi, the QCD 
corrections are smaller when the result is expressed in terms of the running mass. Second, 
the fj, dependence is flatter if the running mass is used. For small values of \i one further 
observes a strong decrease (increase) of the LO (NLO) prediction for cr tot expressed in terms 
of the running mass. These unphysical features occur because both the coupling and the 
running mass become large as /i becomes smallf]. The decrease of the LO result simply 
reflects the fact that of ot goes to zero as z approaches the threshold value z = 1/4. The 
fact that the NLO result in both mass renormalization schemes is rather stable when /i is 
varied as well as the good agreement of the predictions in the two schemes indicates that 
the contributions of order a 2 are small; this is indeed the the explicit calculation of 

shows. 

We now turn back to the discussion of the fraction r(x cut ). The renormalization scale 
dependence of the NLO coefficient B in the pole mass scheme is easily obtained: 

B(x cnt , z on , fi) = B(x cuu z on , y/s) - A(x cnt , z on )Po In (^) , (34) 

4 I thank M. Spira for explaining to me this method of switching between two mass renormalization 
schemes. 

5 For example, at fi = 25 GeV we have m t (/i = 25 GeV) « 199.3 GeV. 



Az(^) 



2C F a s (fj) 



l--ln 

4 



z(n)s 
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with j3 = (llNc — 2n/)/3 and rif = 6. In analogy to (|33|) we may express the result for 
r ( x cut) i n terms of the MS mass at the scale \i by writing 

r(x cut ,z on ,n) = r(x cut ,z(n) + Az(fx),fi) 

= ^A(z^)) + Az(^(z^)) 

+ (^) B( Xeat ,z(n),n) + 0(aS). (35) 

Figs. 10 and 11 show the LO and NLO results for r at x cnt = 0.1 as a function of \i both in 
the pole mass and the running mass scheme at ^fs = 0.5 TeV and y/s = 1 TeV, respectively. 
At \fs = 0.5 TeV, we find a rather large difference of the NLO results in the two schemes at 
scales n ~ \fs. The additional gluon in the final state carries (at x cut = 0.1) a momentum 
of at least 25 GeV, and the allowed phase space depends rather strongly on the top quark 
massQ, which explains the strong /^-dependence of the LO result in terms of the running 
mass. At NLO the //-dependence is drastically reduced, but r still increases with /x. The 
gain in phase space due to the decrease of m t (/i) overcompensates the decrease of a s (fJ>) 
even at NLO. This suggests that the running mass should not be used in the fixed-order 
prediction for the fraction r at y/s = 0.5 TeV and x cut > 0.1. At these values, we are 
too close to the kinematic threshold of the tig final state where the running mass is an 
"unnatural" parameter. 

At \fs = 1 TeV the differences of the NLO results for the two mass definitions are smaller 
as compared to the case of y/s = 0.5 TeV. Also, the NLO results only weakly depend on 
the choice of /t, especially if the running mass is used. The two scales m t and \fs are now 
rather far apart, and the running mass parameter becomes preferable. For even higher cm. 
energies, the fraction r as defined in (|j) eventually suffers from large contributions in the 
collinear region; as z — > 0, the result in nth order is dominated by terms ~ (a s ln(z)) n , 
which should either be resummed or avoided by modifying the definition of r. One could for 
example, after selecting samples of tiX events with 1 — x t i> x cut) use in addition standard 
jet algorithms to define collinear safe cross sections. 

An intriguing question is whether a measurement of the fraction r allows a direct de- 
termination of the value of the running mass parameter of the top quark at high energies. 
(Analogously, NLO results for three-jet fractions involving b quarks |Tj| 0, [U| have been 



used to extract a value for m&(// = mz) from the high-statistics LEP []19| and SLD |20|, |21 
data.) We consider here study the result for r(x cut = 0.1) at yfs = 1 TeV. As was 

discussed above, for these values the fraction r is perturbatively well under control when 
expressed in terms of the running mass. We now vary the value mj(/i = 1 TeV) between 
140 and 160 GeVQ and compute r {xcut — 0.1) at LO and NLO. The results are shown in 
Fig. 12. At NLO, 

K^cut = 0.1) decreases by about 4% when the running top quark mass is 

6 The running top quark mass decreases from 199.3 GeV to 155.3 GeV as fj, is varied between 25 and 500 
GeV. 

7 The value obtained from the renormalization group evolution is m t (fj, = 1 TeV) = 148.6 GeV. 
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changed from 140 to 160 GeV. If a measurement of r will be possible with an error of ±1%, 
the running top quark mass at fi — 1 TeV could be determined up to ±5 GeV. A statistical 
error of 1% on r is realistic with the envisioned high luminosity of a future linear collider 
operating at y/s = 1 TeV, but whether or not the systematic error can be kept that small 
remains to be seen. 

Although the sensitivity of the fraction r on = 1 TeV) seems to be rather poor, a 
direct determination of the running top quark mass at such a high scale would provide a 
nice test of perturbative QCD in the following way: It is expected that a very precise value 
of the top quark mass can be obtained from the threshold scan of a tot . Since the pole mass 
suffers from (nonperturbative) infrared ambiguities, improved mass definitions have been 
proposed in this context, the so-called potential-subtracted [|, || and IS 0] mass. After 
extracting these masses from the threshold scan, they can be converted into the MS mass, 
evolved to /x = 1 TeV, and compared to the direct measurement suggested above. 

5 Conclusions 

At a future high-energy and high-luminosity e + e~ linear collider, the production of top 
quark pairs in association with one or more additional hard partons will be an interesting 
new testing ground for perturbative QCD. For example, a measurement of the fraction r 
of tiX events with respect to all ti events will provide a powerful "flavour independence 
test" of the strong interactions. In this paper we have studied the reaction e + e~ — > tiX 
to order a 2 . The fraction r has been defined as a function of a preset cut parameter x cut 
by demanding 1 — x t i > x cvA where x t i = (h + h) 2 / s is the scaled invariant mass square 
of the ti pair. At LO, where 1 — x t i is equal to the scaled gluon energy x g = 2E g /y/s, we 
have calculated this fraction analytically. The QCD corrections to r(x cut ), which we have 
evaluated for cm. energies y/s = 0.5 TeV and y/s = 1 TeV, are large and positive. We have 
further computed distributions of the differential cross section for e + e~ — > tiX with respect 
to the scaled top quark energy, the variable 1 — x t i, and the cosine of the angle between t 
and i. A measurement of such distributions will be interesting for instance in the context 
of searches for possible anomalous top quark couplings, and the QCD corrections have to 
be included in such analyses. Finally we have studied the dependence of our results on the 
renormalization scale /i and compared the pole mass and the running mass renormalization 
schemes. We have found that at y/s = 0.5 TeV the phase space for tig production depends 
(for a hard gluon) rather strongly on the value of the top quark mass parameter. We have 
concluded that for this (or smaller) cm. energies the fixed order prediction for the fraction 
r should not be expressed in terms of the running mass. At = 1 TeV, fn t (n) becomes the 
preferable mass parameter. This is reflected in the improved stability of the NLO results 
under variations of the scale /i. If a measurement of r at y/s = 1 TeV is possible to an 
accuracy of 1%, the running mass m t (/i = 1 TeV) could be directly determined up to about 
±5 GeV, offering a nice consistency check of the renormalization group evolution of m t (fj,) 
as predicted by perturbative QCD. 

In this paper we have concentrated on studying the fraction of tiX events and their 
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distributions with respect to energies of and angles between final state particles only. Future 
work will include the investigation of observables that involve the orientation of the final 
state with respect to the electron beam, like forward-backward asymmetries or so-called 
"event handedness" correlations [ffO . 
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TABLE CAPTIONS 



Table 1. Leading order and next-to- leading order coefficients of the fraction r(x cut ,/i) as 
denned in Eq. fl27D for two different center-of-mass energies and x cut = 0.1. The renormal- 
ization scale \i is set to \/s, and a top quark pole mass of m t = 175 GeV is used. 

FIGURE CAPTIONS 

Fig. 1. Kinematically allowed region in the (x,x) plane for a tig final state with m t = 175 
GeV. The enveloping light-coloured triangle is the allowed region for massless quarks. Lines 
indicate cuts on the scaled gluon energy. 

Fig. 2. Fraction r(x cut ) (defined in (H)) at leading order in a s for x cut = 0.02,0.1,0.2 and 
m t = 175 GeV as a function of the center-of-mass energy yfs. The running of a s is included. 
Fig. 3. Fraction r(x cut ) at leading (LO) and next-to-leading order (NLO) as a function of 
a^cut at y/s = 0.5 TeV. The renormalization scale /i is set to y/s. 
Fig. 4. Same as Fig. 3, but for y/s = 1 TeV. 

Fig. 5. Distribution 1/ Ototda?,/ dx at LO and NLO for y/s = /i = 0.5 TeV and x cut = 0.1. 
The scaled quark energy is defined in Eq. (|3|), and the cross section cr 3 is defined in Eq. 

Fig. 6. Same as Fig. 5, but for l/a tot d<7 3 /d(l — x t t) with x t i defined in Eq. (|3p. 

Fig. 7. Same as Fig. 5, but for l/a tot da 3 /d(cos9 t i) with cos 6^ defined in Eq. ([T9|). 

Fig. 8. LO and NLO results for the total ti cross section in units of the point cross section 

defined in Eq. (^|) at a cm. energy of yfs = 0.5 TeV. The renormalization scale /i is varied 

between 25 GeV and y/s, and results obtained in the pole mass and in the MS (running) 

mass renormalization schemes are compared. 

Fig. 9. Same as Fig. 8, but for ^/s = 1 TeV. 

Fig. 10. Same as Fig. 8, but for the fraction r at x cut = 0.1. 

Fig. 11. Same as Fig. 10, but for ^fs = 1 TeV. 

Fig. 12. Dependence of the fraction r at x cnt = 0.1 and yfs = \i = 1 TeV on the value of 
the running top quark mass m t at the scale [i = 1 TeV. The dashed line is the LO result, 
and the solid line is the NLO result. 
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